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Abstract 


This  paper  discusses  optimization  techniques  for  the  determination  of  complex  per- 
mittivity and  permeability  in  transmission  lines.  The  traditional  theoretical  model  using 
scattering  parameters  is  extended  into  a mathematical  regression  model  that  can  be  solved 
with  widely  accepted  numerical  techniques.  This  new  model  produces  accurate  primary 
mode  results  for  the  samples  tested  including  nonmagnetic  and  magnetic  materials  with 
high  dielectric  constants.  An  extension  of  the  model  includes  responses  due  to  higher  order 
modes.  The  general  model  determines  parameters  to  specify  the  spectral  functional  form  of 
complex  permittivity  and  permeability  and  is  capable  of  small  corrections  to  independent 
variable  data  including  angular  frequency,  sample  length,  sample  position,  and  cutoff  wave- 
length. The  method  provides  reliable  determination  for  both  low  and  high  permittivity 
materials. 

Key  words:  higher  order  modes,  microwave,  permeability,  permittivity,  primary  mode, 
orthogonal  distance  regression,  scattering  matrix. 
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1.  Introduction 


A constrained  nonlineaLT  optimization  procedure  is  presented  for  the  determination  of  complex 
permittivity  and  permeability  spectra  from  scattering  parameter  (S-parameter)  data  taken  from 
an  automatic  network  analyzer  (ANA).  The  procedure  has  been  used  successfully  for  reliable 
characterization  of  permittivity  and  permeability  of  many  different  test  samples.  In  addition,  it 
provides  a basis  for  the  analysis  of  multi-mode  field  data  and  for  the  determination  of  experi- 
mental systematic  uncertainty. 

Previous  work  in  this  area  involved  the  determination  of  permittivity  and  permeability  on 
a point-by-point  basis  with  explicit  or  implicit  solution  of  a system  of  nonlinear  scattering 
equations  at  each  particular  frequency  (see  [6],  [3]).  Inaccurate  results,  however,  may  arise 
when  numericzd  singularities  occur  at  frequencies  corresponding  to  integer  multiples  of  one  half 
wavelength  of  the  material.  The  same  system  of  nonlinear  scattering  equations  is  used  in  this 
study.  Here  they  are  solved  in  the  sense  of  least  squares  over  the  entire  range  of  measurement  by 
determining  the  best  Laurent  series  approximations  to  permittivity  and  permeability  consistent 
with  linearity  cind  causality  constraints.  Points  of  singularity  may  be  de-emphasized  to  lessen 
the  effect  of  highly  uncertain  data  points. 

This  effort  determines  complex  permittivity  and  permeability  from  two-port  S-parameter 
data  using  the  primary  (or  fundamental)  mode  field  behavior  in  various  materials.  Physical 
measurements  of  the  S-parameter  data  are  made  with  an  automatic  network  analyzer.  Funda- 
mental mode  S-parameter  relationships  are  then  used  to  solve  for  permeability  and  permittivity, 
and  e(a;)  respectively,  as  a function  of  the  angular  frequency  w. 

The  general  functional  form  for  p and  e was  decided  upon  after  evaluating  several  different 
polynomial  and  trigonometric  relations.  The  evaluation  criteria  for  the  functional  form  for 
permittivity  and  permeability  was  based  on  the  reproducibility  of  the  S-parameter  data  in  terms 
of  the  minimum  total  least  square  approximation  to  these  data.  The  best  overall  approximation 
was  selected.  The  general  functional  form  for  fj,  and  e with  acceptable  results  involves  only  the 
first  two  terms  of  a Laurent  series. 

The  optimization  approach  in  this  research  is  an  implicit  function  regression  model.  This 
model  is  solved  by  the  orthogonal  distance  regression  package  ODRPACK  [2].  This  estimation 
package  allows  for  adjustments  in  input  parameters  to  compensate  for  measurement  uncertain- 
ties. Adjustments  here  are  limited  to  the  sample  length,  scimple  position  in  the  waveguide 
measurement  fixture,  and  the  cutoff  wavelength.  Orthogonal  distance  regression  is  intended  to 
compensate  for  slight  uncertainties  in  the  independent  variable  (angular  frequency,  a;)  as  well 
cis  the  dependent  variable  (observed  S-parameter  data). 

The  physical  model  is  outlined  in  §2,  with  the  related  mathematical  model  discussed  in  §3. 
Numerical  considerations  tire  covered  in  §4  and  conclusions  and  future  directions  discussed  in 

§5. 
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2.  Scattering  Parameter  Relations 


The  equations  described  below  relate  the  measured  two-port  scattering  parameters  (S- 
parameters)  to  the  permittivity  and  permeability  of  the  material.  First,  in  order  to  develop 
the  scattering  equations,  the  following  notation  is  used.  Let 


^ — ^*R^o 

and 

= Wr  - = fJ'*RfJ'o, 

the  permittivity  and  permeability  of  a sample  material,  where  Cq  and  fio  are  the  permittivity  and 
permeability  of  a vacuum,  and  and  are  the  relative  complex  permittivity  and  permeabihty. 
Next  let  Cyac  and  ciab  be  the  speed  of  light  in  a vacuum  and  the  laboratory,  respectively,  and 
for  a given  frequency  /,  let 

0}  = 2k  f, 

be  the  corresponding  angular  frequency.  Then 


7 = 


and 


Vt?-©’ 


represent  the  propagation  constants  in  the  material  and  air,  respectively,  where  j = \/— 1 and 
Aci  is  the  cutoff  wavelength  in  the  waveguide  measurement  fixture,  where  the  subscript  1 refers 
to  the  fundamental  mode.  The  expression  for  the  transmission  coefficient  z is 

z = exp{—'iL), 

where  L is  the  sample  length.  The  reflection  coefficient  is 


r = 


or 


r = 


£j 
Clab 
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Figure  2.1:  Dielectric  Sample  in  Waveguide 


It  is  assumed  that  the  total  length  of  the  sample  holder  is 

Lair  = L + Li  -{■  1-2,  (2.1) 

as  shown  in  figure  2.1,  where  Li  and  L2  are  the  distances  from  the  calibration  reference  planes 
to  the  sample  faces  for  ports  1 and  2,  respectively. 

For  a two-port  device  the  expressions  for  the  measured  S-parameters  are  obtained  by  the 
solution  of  a related  boundary  value  problem.  The  explicit  expressions  for  the  scattering  relations 
of  the  fundamental  mode  are  assumed  to  functions  of  and  axe  given  by 


5ii(A,J 

= 

■S'l2(-^ci) 

= RiR 

'52l(Aci) 

= RiR 

•522('^Ci) 

= Rl 

nini!)' 

1 - r2z2 

^(i-r^) 


1 - r2z2 
^(i-r^) 


1 - r202 


1 - r2z2 


(2.2) 

(2.3) 

(2.4) 

(2.5) 

(2.6) 


where  Ri  — exp(— 70^1)  and  R2  = exp(— 70^2)  are  the  reference  plane  transformation  expres- 
sions. 
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3.  The  Mathematical  Model 


3.1.  e^(w),  and  the  Regression  Model 


The  mathematical  problem  can  be  stated  as  that  of  finding  parameters  to  a prespecified  func- 
tional form  for  the  complex  functions  and  e'^{u>).  S-parameter  data  acquired  from  the 

ANA  for  selected  samples  of  materials  provide  complex  values  for  each  of  5ii,  S21,  S12,  and  S22, 
at  n different  frequencies  ranging  from  1 to  18  GHz.  The  model  determines  the  parameterization 
of  and  that  best  reproduces  simultaneously  the  four  S-parameters  for  the  n observations, 
given  the  sample  data,  two  reference  plane  positions,  and  the  sample  length. 

The  general  form  for  /x^(a;)  and  uses  terms  from  the  Laurent  series: 


fH 


(1  -t-  biu>y  ’ 


(3.1) 


where  a,-  and  bi  are  complex  scalars  in  the  term.  The  solution  procedure  here  is  an  implicit 
function  regression  and  corresponds  to  the  parameterization  of  two  terms  in  a truncated  Laurent 
series.  Note  that  eq  (3.1)  automatically  satisfies  Kramers- Kronig  relations  (see,  for  example,  [l] 
or  [4])  for  dispersion.  The  functional  form  for  /x^(a;)  and  c^(w)  follows: 


= f(<^) 
4(g^)  = f(^) 


, ^2 

1 -f-  BiU}  (1  -h  82^0)^’ 

^2 

1 -)-  Bioj  (1  + B2(^)'^ 


For  notational  convenience,  the  terms  of  the  second  truncated  Laurent  series  Ai,  A2,  Bi,  and 
B2  will  be  referred  to  as  A3,  A^,  B3,  and  B^. 

The  solution  procedure  to  determine  the  complex  parameters  Ai,Bi  is  the  minimization  of 
the  sum  of  the  squares  of  the  uncertainties  between  the  predicted  and  observed  S-parameters, 


min  El’Sn-'f’nl'  + l® 


rj-ip  + isji-p'p  + is. 


(3.2) 


Kk-1 


where,  for  i,j  = 1,2,  5,^-  is  the  observed  5,y  scattering  parameter  at  frequency  oJk  and 


is  the  corresponding  predicted  scattering  parameter.  We  let  l^j  = -f  ^(z)^  represent 

the  absolute  value  of  a complex  scalar  z.  Note  that  eq  (3.2)  is  equivalent  to  the  minimization 
of  the  sum  of  the  squares  of  the  real  and  imaginary  parts  of  each  of  the  S-parameters,  i.e., 

i. 


min  I E E E - Pf’if  + " P^i)'' 

k=l i=l 


(3.3) 
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The  model  reported  in  §4.4  uses  the  formulation  in  eq  (3.3). 


S.2.  Adjustments  to  the  Model  Inputs 

In  standard  ordinary  least  squares  regression  models,  the  observed  responses  are  assumed  to 
contain  some  uncertainties  either  produced  by  the  phenomena  under  examination  or  introduced 
by  the  device  that  measures  the  events.  In  addition,  certain  independent  and  dependent  variable 
pairs  may  not  be  as  reliable  as  others  due  to  an  increased  variance  in  the  uncertainties  in  the 
dependent  variable  for  particular  values  of  the  independent  variable.  The  application  of  a 
reduced  weight  to  points  of  questionable  reliability  provides  the  modeler  with  a means  to  de- 
emphasize  such  points  to  find  a more  appropriate  regression  solution. 

Usually  the  independent  variable  can  be  controlled,  and  the  precise  value  of  each  of  the 
observations  is  well  known.  An  orthogonal  distance  regression  model  provides  the  modeler  with 
the  additional  ability  to  assume  that  the  independent  variable,  in  this  case  frequency,  may 
contain  some  uncertainty  eis  well.  Allowances  for  this  type  of  uncertainty  can,  in  some  cases, 
greatly  improve  the  approximation.  For  this  particular  model  and  the  samples  tested  in  this 
study,  the  uncertainty  in  the  independent  variables  is  sufficiently  small  to  allow  the  modeler  to 
cissume  that  an  ordinary  least  squares  approximation  provides  an  adequate  solution. 

Other  model  parameters  such  as  sample  length,  sample  position  in  the  waveguide,  and  cutoff 
wavelength  are  sufficiently  sensitive  to  require  slight  perturbation.  Each  of  these  inputs  is 
required  by  the  S-parameter  equations  and  is  considered  to  be  known  by  the  modeler.  Since 
these  inputs  eire  not  always  known  exactly,  each  may  be  perturbed  slightly,  as  determined  by 
the  problem  solver,  to  improve  approximation.  This  allows  the  user  to  adjust  for  measurement 
uncertainty. 

In  particular,  since  incorrect  specification  of  the  sample  position,  Li,  in  waveguide  affects 
the  value  of  phase  in  the  reflected  S-parameter  data,  an  additional  parameter,  is  included 
in  Ri, 

Ri  = exp(-7o[Xi  -F  I3li]). 

Similarly,  uncertainties  in  Lair  ^d  L2  are  represented  by  analogous  means  with  parameters 
and  respectively.  With  eq  2.1,  the  total  length  L of  the  sample  is  completely  determined 
by 

L = {Lair  + 0La)  ~ (-^'1  + Pli  + L2  + Pl^) 

and  is  pcirameterized  by  the  values  of  P^a,  Pbiy  and  Pl2- 

Another  parameterized  correction,  Px,  is  included  in  the  cutoff  wavelength  Aci  as  + Px- 
Inaccuracies  in  the  waveguide  dimensions  due  to  the  milling  process  can  affect  the  value  for  the 
cutoflF  wavelength.  In  addition,  the  higher  order  model  discussed  below  requires  an  additional 
wavelength  cutoff  value,  Ac,,  for  a higher  mode  solution.  Since  it  is  not  known  a priori  which 
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higher  modes  will  be  present,  a variation  in  Acj  from  0 < , enables  the  solution 

procedure  to  find  the  value  for  Acj  that  best  improves  the  S-parameter  approximation  for  the 
higher  order  terms. 

One  can  formulate  an  invariant  model  with  respect  to  reference  planes  for  the  problem 
discussed  here,  as  suggested  in  [l],  to  remove  uncertainty  in  Li  and  in  L2.  This  approach  uses 
eq  2.3  or  2.4,  and  the  determinant  of  the  S-matrix, 

det{S)  = 5ii5'22  — -S'i2<S'2i, 

and  solves  for  the  parameters  Ai, . . . , A^^Bi . . , B^.  Simplification  of  det[S)  yields  a formula 
with  neither  Li  nor  L2-  This  reduced  model  remains  dependent  on  a precise  knowledge  of 
both  Lair  and  L and  the  cutoff  wavelength.  Ad-  This  approach  can  be  used  interchangeably 
with  the  original  model  that  contains  Li  and  L2-  The  original  model  uses  twice  the  number  of 
observations  and  seems  to  produce  more  accurate  approximations  to  the  S-parameter  data. 

S.3.  Higher  Order  Modes 

In  samples  with  a high  dielectric  constant,  the  observed  S-parameter  data  may  exhibit  responses 
due  to  modes  other  than  the  fundamental  mode.  These  higher  modes  are  the  result  of  resonances 
of  the  fundamental  mode  in  the  material.  Earlier  work  does  not  include  this  information  in 
the  computation  of  and  as  higher  modes  are  ignored.  In  this  section  we  describe  an 
enhancement  to  the  previously  defined  model  that  does  include  the  higher  order  response  data. 

Because  of  the  similarity  between  the  response  of  the  primary  mode  and  the  higher  modes 
that  is  visually  apparent  (see  §4.4),  an  additional  term  is  added  to  each  of  the  four  equations  (see 
2. 2-2. 5).  The  new  numerical  model  includes  higher  order  mode  structure  in  the  S-parameter 
approximation  and  is  a simple  extension  of  the  primary  mode  model  described  in  §2.^  The 
terms  that  approximate  the  higher  order  modes  are  identical  to  the  primary  mode  term  with 
the  exception  that  each  term  is  scaled.  For  the  evaluation  of  the  higher  order  terms,  all  input 
parameters  are  unchanged  except  for  the  parameter  for  the  higher  mode  cutoff  wavelength,  Ac,, 
that  is  allowed  to  decrease. 

The  explicit  forms  of  the  new  set  of  equations  use  the  eqs  (2. 2-2. 5),  and  are  denoted 
as  Dll,  Di2,  D21,  3^nd  Z?22-  Associated  with  each  higher  order  term  is  a scaling  parameter, 
Pi,. . . , P4.  In  particular,  to  include  the  primary  mode  and  one  higher-order  mode  the  following 
approximate  model  is  used: 

Dll  = ‘5ii('^ci)  + ^i5ii(Ac3),  (3.4) 

D21  — •S2i(Aci)  + /3252i(Ac3),  (3.5) 

*The  justification  of  new  terms  added  to  each  of  (2. 2-2. 5)  is  based  solely  on  empirical  evidence  found  when 
solving  for  the  primary  mode. 
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5i2(Aci)  + PzSi2{^ct)t 
‘S'22(-^Ci)  + pAS^-zi^Ci)^ 


(3.6) 

(3.7) 


£>12  — 

Z?22  = 

where  Ac,  is  determined  by  the  problem  solver. 

3.4.  The  Initial  Solution  Procedure 

In  the  Nicolson-Ross-Weir  procedure  ([5],  [7])  the  equations  for  the  scattering  parameters  are 
combined  to  allow  the  system  of  equations  to  decouple.  This  decoupling  yields  an  explicit 
equation  for  the  permittivity  and  permeability  as  a function  of  the  scattering  peirameters  on  a 
point-by-point  basis.  This  solution  procedure  is  the  basis  of  the  commonly  used  techniques  for 
obtaining  permittivity  and  permeability.  Unfortunately,  these  equations  are  not  well  behaved 
for  low-loss  materials  at  frequencies  that  correspond  to  integer  multiples  of  one  half  wavelength 
in  the  sample. 

The  Nicolson-Ross-Weir  procedure  provides  a good  initial  approximation  to  and 

B^.  All  other  parameters  are  initialized  to  zero.  The  estimated  values  for  permittivity  and  per- 
meability are  determined  on  a point-by-point  basis  by  frequency.  The  corresponding  scattering 
parameters  are  computed  with  these  values  and  then  compared  to  the  observed  values.  The 
computed  values  for  and  that  provide  the  closest  agreement  between  the  observed  and 
predicted  S-parameter  data  are  used  as  the  initial  values  to  the  regression  model. 
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4.  Numerical  Considerations 


4.1.  ODRPACK  Orthogonal  Distance  Regression  Package 

Briefly,  the  solution  procedure  ODRPACK  [2]  is  trust  region  Levenberg-Marquaxdt  approach. 
This  type  of  trust  region  approach  adaptively  determines  the  region  in  which  the  linear  approx- 
imation closely  resembles  the  nonlinear  model.  The  procedure  allows  both  the  ordinary  least 
squares  model,  in  which  the  uncertainties  are  assumed  to  be  only  in  the  dependent  variable, 
and,  the  orthogonal  distance  regression  model,  where  uncertainties  also  exist  for  the  indepen- 
dant variables. 

First  order  derivatives  for  the  Jacobian  matrices  can  be  numerically  approximated  (finite 
difference  approximation),  or  can  be  user-supplied  analytical  derivatives.  The  procedure  per- 
forms automatic  scaling  of  the  variables  if  necessary,  as  well  as  determining  the  precision  of  the 
model  in  terms  of  machine  precision.  The  ODRPACK  includes  many  other  features  that  assist 
the  user  in  the  modeling  process.  The  model  automatically  determines  the  number  of  usable 
digits,  checks  analytical  derivatives  provided  by  the  user,  and  automatically  selects  many  of  the 
input  parameters  for  the  user,  if  desired. 

Iterations  are  stopped  by  ODRPACK  when  any  one  of  three  stopping  criteria  are  met.  Two 
of  these  indicate  that  the  iterations  have  converged  to  a solution.  Sum-of-squares  convergence 
indicates  that  the  change  in  sum-of-squares  observational  uncertainty  is  sufficiently  small,  and 
parameter  convergence  indicates  that  the  change  in  the  estimated  parameters  is  sufficiently 
small.  The  third  stopping  criterion  is  a limit  on  the  number  of  iterations. 

4.2.  Initial  Conditions 

Many  of  the  input  options  for  the  problem  solver  can  be  set  to  their  default  values,  as  is  done  in 
this  model.  The  most  significant  input  parameters  for  modeling  permittivity  and  permeability 
are  the  initial  values  for  Ai,  A3,  Bi  and  B3.  Sensitivity  to  the  initial  solution  for  these  parameters 
is  discussed  below  and  the  selection  of  initial  settings  of  the  starting  values  is  covered  in  §3.4. 
When  higher  modes  are  included,  the  solution  for  the  primary  mode  is  used  cis  the  initial  guess 
for  the  higher  mode  model.  For  all  of  the  parameters  that  define  the  physical  model  except 
those  for  and  mentioned  above,  the  standard  laboratory  values  are  used.  All  additional 
parameters  are  initialized  to  zero. 

4.3.  Sample  Characteristics 

A total  of  seven  samples  were  modeled  to  determine  permittivity  and  permeability  with  ANA 
two-port  S-par£imeter  data.  The  sample  characteristics  appear  in  table  4.1.  Figure  4.1  contains 
selected  plots  of  the  observed  ANA  S-parameter  data  versus  frequency  for  sample  #6. 


8 


Table  4.1:  Sample  chciracteristics 


Material 

Name 

Identifier 

Length 
in  meters 

Li 

L2 

Initial 

Initial 

l^*R 

1 

cross-linked  polystyrene 

rexa240889 

2.407x10-2 

0 

1.3208x10-3 

(2.53,0.002) 

(1.0) 

2 

1723  glass 

172a240889 

1.015x10-2 

0 

1.5250x10-2 

(6.15,0.04) 

(1.0) 

3 

loaded  polymer 

112a050290 

2.540x10-2 

0 

0.0000000 

(5.75,0.23) 

(1.6,0.1) 

4 

yttrium  iron  garnet 

YIG 

1.766x10-2 

0 

7.7380x10-3 

(10,0.2) 

(1,0.4) 

5 

nickel  ferrite 

ttlal20490 

1.013x10-2 

0 

1.5263x10-2 

(11.5,0.1) 

(0.85,0.1) 

6 

barium  titanate  mix#l 

barium 

7.632x10-3 

0 

1.7768x10-2 

(265,1) 

(1.0) 

7 

barium  titanate  mix^2 

french 

2.427x10-2 

0 

1.1252x10-3 

(105,1) 

(1.0) 

4.4.  Numerical  Results 

The  first  results  reported  are  those  for  the  primary  mode  model  for  all  seven  samples  described 
in  §4.3.  The  first  set  of  plots  (figure  4.2)  include  the  observed  S-parameter  data  (dots)  for  the 
YIG  sample  (scimple  #4)  from  the  ANA  overlaid  by  the  predicted  data  (line)  found  by  the 
model.  The  corresponding  residual  plots  ^ for  this  sample  appear  in  figure  4.3. 

For  the  first  four  samples  the  predicted  and  observed  data  are  nearly  identical.  The  residual 
plot  for  cross-linked  polystyrene  (sample  #1)  as  shown  in  figure  4.4  reveab  only  the  systematic 
uncertainty  due  to  the  ANA. 

Figure  4.5  illustrates  a S-parameter  primary  mode  solution  for  the  bcirium  titanate  with 
a high  dielectric  constant  ^lnd  its  corresponding  residual  plot.  For  three  of  the  samples  with  high 
dielectric  constants,  the  new  model  produces  responses  for  the  primary  and  one  of  the  higher 
modes.  See  figure  4.6  for  the  S-parameter  real  and  imaginary  components  of  ^21  for  the  sample 
exhibited  in  figure  4.5  with  the  model  for  higher  order  modes. 

For  samples  ^^5,  #6,  and  #7  the  predicted  S-parameter  data  provide  realistic  primary  mode 
responses,"  although  the  residual  plots  for  these  samples  indicate  that  higher  modes  are  present. 
With  the  higher  order  model  described  in  §3.3  and  the  new  solution  found,  the  problem  is  re- 
solved for  the  last  three  samples.  The  new  solution  specifies  the  additional  parameters  Pi,. . . , 
Figure  4.7  shows  the  results  from  the  higher  order  model  for  sample  #5.  Additional  work  in  this 
area  suggests  that  this  model  is  only  an  approximate  theoretical  model  for  the  higher  modes. 

“residual  plot”  denotes  a plot  of  the  difference  between  the  predicted  and  observed  S-parameter  data  over 
the  range  of  observations. 
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Figure  4.1:  Barium  Titanate  #1  (a)  (b)  9[<5ii],  (c)  9i[52i],  (d)  9f[52i] 

4.5.  ANA  Systematic  Uncertainty 

The  difference  of  the  predicted  S-parameter  and  the  observed  values  for  samples  with  low  di- 
electric values  revealed  systematic  uncertainty.  (For  the  high  dielectric  s^ples,  other  sources 
of  uncertainty  including  higher  mode  responses  dominate  the  ANA-induced  uncertainty.)  Ad- 
ditional tests  revealed  that  the  source  of  the  uncertainty  is  not  related  to  the  material  tested  in 
the  waveguide.  In  fact,  uncertainties  in  the  S-parameter  data  for  the  cross-linked  polystyrene 
sample  closely  resemble  the  S-parameter  data  for  an  empty  waveguide.  Attempts  to  further 
identify  and  remove  the  systematic  uncertainty  are  underway.  Figure  4.4  contains  the  data 
for  an  empty  waveguide  and  also  the  residual  plot  for  the  cross-linked  polystyrene  sample. 

4.6.  Permittivity  and  Permeability  Estimates 

The  estimates  of  permeability  and  permittivity  for  the  Vcirious  samples  are  defined  by  the  pa- 
rameters Ai, . . . A4,  Bi,. . .B4  as  a function  of  frequency.  Permittivity  values  of  several  of  the 
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Figure  4.2:  Yttrium  iron  garnet  predicted  and  observed:  (a)  9i[5i2])  (b)  9f[5i2] 

samples  are  equal  to  (1,0)  (see  table  4.1).  For  these  samples  the  permittivity  was  allowed  to 
vary,  with  the  problem  solver  required  to  determine  the  (1,0)  value.  Slight  veiriations  in  the 
value  of  permittivity  are  apparent,  although  the  deviation  from  the  true  value  is  small.  The 
solution  procedure  also  provides  the  standard  deviation  and  confidence  intervals  for  each  of  the 
estimated  parameters. 

4.7.  Robustness  of  the  ODRPACK  Procedure 

The  robustness  of  a mathematical  procedure  is  its  ability  to  find  a locally  optimal  solution 
from  a varied  initial  condition.  The  robustness  of  the  entire  permeability  and  permittivity 
procediue  depends  on  the  robustness  of  the  ODRPACK  procedure  and,  more  significantly,  the 
robustness  of  the  mathematical  model.  The  existence  of  alternative  optima  in  the  mathematical 
model  can  limit  the  range  of  the  initial  conditions  to  produce  a particular  solution.  In  addition, 
singularities  in  the  absence  of  alternative  optima  may  force  the  solution  procedure  to  fail  to 
determine  directions  of  improvement  and  cause  a premature  termination. 

For  the  samples  in  this  study,  variability  in  the  robustness  of  the  procedure  depended  on 
the  sample  studied.  For  the  materials  with  small  dielectric  constant,  the  procedure  readily 
determined  the  correct  solution  for  a variety  of  initial  conditions.  For  materials  with  higher 
dielectric  constant,  the  procedure  often  converged  quickly,  although  the  existence  of  alternative 
optima  in  the  mathematical  model  often  required  the  repetition  of  tests  with  varied  initial 
conditions  before  an  acceptable  solution  was  found. 

In  paurticular,  for  the  cross-linked  polystyrene  sample,  initial  values  for  and  were 
set  to  [(1, 0),  (1, 0)]  and  were  constant  over  the  entire  frequency  range.  The  correct  solution 
((1,0),  (2.517,0.0018)1  was  found  after  50  iterations  of  the  solution  procedure.  Initial  conditions 
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Figure  4.3:  Yttrium  iron  garnet  residual  plots:  (a)  9i[<Si2],  (b)  9l[5i2] 

of  [(1, 0),  (4.0, 0.01)]  failed  to  determine  the  correct  solution,  while  altering  the  imaginary  part 
of  to  1 X 10“'*  resulted  in  the  correct  final  solution  once  again.  For  material  with  high 
dielectric  constants,  the  range  was  smaller  in  percentage  but  Icirger  in  amount.  For  example,  for 
the  bcirium  titanate  #1,  initial  conditions  of  [(1.0, 0.0),  (200.0, 1.0)]  produced  the  correct  value 
of  [(1,0),  (269.0, 1.70)]  in  59  iterations. 

One  should  note  that  the  solution  procedure  can  change  the  value  of  the  estimated  parameters 
to  the  model  by  large  amounts  in  the  early  iterations.  Hence  it  is  possible  that  for  cross-linked 
polystyrene  an  initial  solution  farther  from  the  desired  solution  may  in  fact  find  the  correct 
solution  merely  because  of  the  solution  steps  taken  by  the  procedure. 
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Figure  4.4:  (a)  Residual  Plots  of  9[522]  for  cross-linked  polystyrene  (b)  9(522]  for  an  empty 
waveguide 
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Figure  4.5:  Barium  titanate  #1  (a)  Predicted  and  Observed  8?(52i]  (b)  residual  plot 
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Figure  4.6:  Barium  titanate  #1  (a)  Predicted  and  Observed  5i[52i]  (b)  Predicted  and  Observed 
S[52i] 


Figure  4.7:  Nickel  ferrite  (a)  Predicted  and  Observed  3i[5i2]  (b)  Predicted  and  Observed  9[5i2] 
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Figure  4.8:  Selected  plots  of  for  various  samples 
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5.  Conclusions  and  Future  Directions 


The  nonlinear  optimization  procedure  using  two-port  scattering  parameters  determined  permit- 
tivity and  permeability  for  a number  of  samples  with  a high  degree  of  accuracy.  The  added 
capability  that  permits  variations  in  certain  input  parameters  provides  a mechanism  to  adjust 
for  measurement  uncertainties.  The  extension  of  the  model  to  include  higher  order  mode  re- 
sponses is  quite  useful  for  high  dielectric  constant  materiab.  This  procedure  can  be  extended 
to  include  additional  higher  order  mode  responses  by  the  inclusion  of  more  terms  as  described 
in  §3. 
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